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ABSTRACT 

The recent measurement of the Shapiro delay in the radio pulsar PSR J1614— 2230 yielded a mass 
of 1.97±0.04 Mq, making it the most massive pulsar known to date. Its mass is high enough that, 
even without an accompanying measurement of the stellar radius, it has a strong impact on our 
understanding of nuclear matter, gamma-ray bursts, and the generation of gravitational waves from 
coalescing neutron stars. This single high mass value indicates that a transition to quark matter 
in neutron-star cores can occur at densities comparable to the nuclear saturation density only if the 
quarks are strongly interacting and are color superconducting. We further show that a high maximum 
neutron-star mass is required if short duration gamma-ray bursts are powered by coalescing neutron 
stars and, therefore, this mechanism becomes viable in the light of the recent measurement. Finally, 
we argue that the low- frequency (< 500 Hz) gravitational waves emitted during the final stages of 
neutron-star coalescence encode the properties of the equation of state because neutron stars consistent 
with this measurement cannot be centrally condensed. This will facilitate the measurement of the 
neutron star equation of state with Advanced LIGO/Virgo. 

Subject headings: pulsars — neutron star physics 



1. INTRODUCTION 

Neutron stars are associated with the most diverse 
and energetic phenomena in the Universe, from gamma- 
ray bursts to the emission of gravitational waves and 
from periodic millisecond radio signals to month-long 
X-ray outbursts. The strength and even the occur- 
rence of some of these phenomena depend on the 
neutron star mass. For decades, precise dynami- 
cal measurements of the masses of neutron stars re- 
sulted i n highly clustered values aroun d 1.25-1.4 solar 
masses (jThorsett fc Cha krabartvl ll999f ). This paradigm 
has recently changed with the measurement of the neu- 
tron star mass for PSR J1614-2230. 

PSR J1614-2230 is a 3.1 ms radi o pulsar in a 8.7 da y 
orbit around a massive white dwarf (jHessels et al.ll2005[ ). 
Its very high inclination, at 89.17°, allowed the detection 
of a strong Shapiro delay signature in highly accurate 
timing observations with the Green Bank Telescope and 
a precise measurement of the neutron star mass. The 
inferred value of 1.97±0.04 M & is by far the highest ob- 
serve d from any neutron star to date (|Demorest et al.l 
l2010h . 

The high mass of PSR J1614— 2230 provides a lower 
limit on the maximum mass of neutron stars and sets 
the dividing line between neutron stars and black holes 
to at least this value. This has significant implications for 
the fraction of neutron stars that will collapse into black 
holes as a result of mass accretion in X-ray binaries. For a 
maximum neutron star mass that is > 2 Mq , the relative 
frequency of black hole binaries that are formed through 
this evolutionary channel is at most 25% (Pfahl, Rappa- 
port, & Podsiadlowski 2003). Such a small fraction would 
help explain the lack of low-mass black holes in binary 
systems with dynamical mass measurements (jOzel et al.1 



l2010bD . 

The highest neutron star mass that can be supported 
against collapse depends very sensitively on the underly- 
ing equation of state. In particular, if quark, hyperon, or 
boson degrees of freedom are excited at high densities, 
the equation of state softens and cannot support massive 
neutron stars (see, e.g., Lattimer & Prakash 2001). The 
discovery of even a single very high mass neutron star 
can, therefore, provide strong constraints on the funda- 
mental properties of ultradense matter. In this Letter, 
we explore in detail the constraints imposed by the re- 
cent mass measurement on the radii of neutron stars and 
on the properties of quark matter in their interiors. 

The maximum neutron star mass also plays a crucial 
role in determining the outcome of the coalescence of two 
neutron stars. These astrophysical events are thought 
to be the central engines of short-duration gamma-ray 
bursts (see reviews by Nakar 2007 and Lee & Ramirez- 
Ruiz 2007). Furthermore, they are expected to be the 
primary sources of gravitational waves that will be de- 
tected by ground-based observatories such as LIGO, 
Virgo, and GEO600. In §3, we investigate the viability 
of the coalescing neutron-star scenario for short-duration 
gamma-ray bursts in view of the above lower limit on 
the maximum mass of neutron stars. Finally, in §4, we 
discuss the prospect of measuring the equation of state 
of neutron-star matter using ground-based gravitational 
wave observatories. 

2. THE STRUCTURE AND COMPOSITION OF NEUTRON 

STARS 

The high mass of the PSR J1614— 2230, in conjuction 
with the requirement of causality for the neutron star 
matter equation of state, places a strong lower limit on 
its radius, independent of the composition of its interior. 
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Fig. 1. — The maximum neutron star mass as a function of two 
parameters of quark matter when the density at which the transi- 
tion from nucleonic to quark matter occurs is equal to 1.5 times the 
nuclear saturation density. The measurement of a pulsar mass of 
> 1.93 Mq from Shapiro delay observations indicates that, if the 
transition to quark matter occurs at densities that are relevant to 
neutron star interiors, such a massive star can be supported against 
collapse only if the quarks are strongly interacting (04 < 0.63). 



Assuming an equation of state that arises from normal 
interactions in the outermost neutron star layers (where 
the density p is less than a characteristic value po) and 
requiring that the sound speed c s everywhere in the neu- 
tron star core is less than the speed of light c, such that 
dP/dp < c 2 , leads to a firm lower limit on th e radius R 
of a neutron star as a function of its mass M ()Lin dblom 
fl98i IGlendennindl200a IKoranda et all 119971) . This is, 
of course, true only for neutron stars, which are gravi- 
tationally bound and have densities that drop smoothly 
to zero in their atmospheres; the causality limit does not 
apply to strange stars. 

The BPS equation of state (Baym, Pethick, & Suther- 
land 1971) is a standard choice for low-density mat- 
ter. Assuming it is accurate up to a density po = 
3 x 10 14 g cm -3 , which is comparable to the nuclear sat- 

2.7xl0 14 gcm- 3 , 



uration density of atomic nuclei p sa t 
the bound imposed by causality is 



R > 2.83 — =- = 8.3 



M 



1.97M P 



km. 



(1) 



All realistic equations of state for gravitationally bound 
stars predict radii that are either very weakly de- 
pende nt on or decreasing with increasing neutron-star 
mass (Lattimcr & Prakash 2001). The causality bound, 
in conjunction with the recent spectroscopic measure- 
ments of neutron-star radii in X-ray binaries that place 
them at values < 12 km (Ozel, Baym, & Giiver 2010, 
Steiner, Lattimer, & Brown 2010), lead to the very nar- 
row range 8.3 km < R < 12 km for all neutron stars. 

Stronger constraints on the composition of the neutron 
star interior can be placed by requiring that the equation 
of state support a maximum neutron star mass that is 
at least as high as the measured value. In general, the 
appearance of new degrees of freedom at and above the 



nuclear saturation density, such as quarks, hyperons, or 
bosons, softens the equation of state and lowers the max- 
imum mass that can be supported against collapse. 

In Figure 1, we show how the mass measurement for 
PSR J1614— 2230 constrains the quark matter equation 
of state. We use a specific model of quark matter, the 
phenomenological equation of state proposed by Alford 
et al. (2005), which has three parameters: 04, a 2 , and 
the bag coefficient. The quartic coefficient 04 measures 
the degree of interaction between quarks, with 04 = 1 
corresponding to free quark plasma. The quadratic co- 
efficient a 2 depends on the mass of the strange quark 
and the color superconducting energy gap; it is large in 
unpaired quark matter with a heavy strange quark, and 
small in color superconducting (for example color-flavor 
locked) quark matter with a light strange quark. We cal- 
culate the maximum neutron star mass M max for a range 
of values of 0,4 and 02 , while fixing the bag coefficient so 
that the transition from nuclear matter to quark matter 
occurs when the baryon density n of the nuclear matter 
is 1.5n sat , where nuclear density is n sa t = 0.16fm~ 3 . For 
the nuclear matter equation of state, we use BPS (Baym, 
Pethick, & Sutherland 1971) at low densities and APR 
(Akmal, Pandharipande, & Ravenhall 1998), which we 
assume is valid up to a density 2 n sat . We chose the APR 
equation of state because it describes well the nucleon- 
nucleon scattering data and is believed to be accurate 
up to the nuclear saturation density. The excluded re- 
gion in the lower right part of the figure is the region 
where, after the transition from nuclear to quark matter 
at n = 1.5rt sa t, there is another transition back to nuclear 
matter at 1.5n sat < n < 2n sat . Such a transition seems 
unlikely to be physical, so we do not calculate masses in 
this region. 

We see from Figure 1 that values of M max compatible 
with the measurement for PSR J1614— 2230 can only be 
achieved if a 4 < 0.63 and a 2 < 10 4 MeV 2 . In Alford 
et al. (2005), such values were associated with strong 
interactions between the quarks (pushing 04 well below 
unity) and color superconductivity (which allows a 2 to 
remain small by cancelling the contribution to a 2 from 
the strange quark mass). 

3. THE MAXIMUM MASS OF NEUTRON STARS AND 
SHORT-DURATION GAMMA-RAY BURSTS 

Short gamma-ray bursts have characteristic timescales 
< 2 s, which are much longer than the dynamical 
timcscalc for the coalescence of two neutron stars. In 
the coalescing neutron star scenario, the re quired l onger 
timescales can be achieved in two ways (iNakarl 1200 
Lee fc Ramirez-Ruiz! I2007t iDuezl 120101 : iRezzolla et al 
2010f ). The initial remnant may be a differentially rotat 
ing massive neutron star that is temporarily supported 
by centrifugal forces. The dissipation of the differential 
rotation on timescales ~ 1 s leads to the delayed col- 
lapse of the remnant into a black hole and produces the 
gamma-ray burst. In order for this channel to occur, the 
sum of the masses in the initial system has to be larger 
than the maximum neutron star mass for solid body ro- 
tation but less than the maximum mass that can be sup- 
ported by differential rotation. Numerical simulations of 
spinning neutron stars place these limits at ~ 1.2 M max 
and 1.4 M max , respectively, where M max is the max- 
imum mass of a non-spinning neutron star (|Duezl [20101 : 
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Fig. 2. — Allowed regions of the parameter space of the masses of two coalescing neutron stars that leads to a short-duration gamma-ray 
burst. For simplicity, the second neutron star is assumed to be larger than the first neutron star. The left panel shows the result when 
the maximum mass of a non-spinning neutron star is equal to M max = 1.8 Mq and the right panel shows the result for a maximum 
mass of 2.1 Mq. In the right panel, the blue line corresponds to a total mass of Mi + M2 = 1.2 M max , below which a neutron star 
rotating as a solid body can be supported against collapse by centrifugal forces. In both panels, the red line corresponds to a total mass 
of Mi + M2 = 1.4 M max , above which a neutron star cannot be supported against collapse, even if it is rotating differentially. The green 
line corresponds to a mass ratio of M2/M1 = 4/3, above which the initial outcome of the collision is a black hole surrounded by a massive 
torus. Short duration gamma ray bursts can be generated when the outcome of the collision is either the delayed collapse of a supermassive 
neutron star into a black hole, or the prompt collapse of the two stars into a black hole surrounded by a massive torus. The allowed region 
of the parameter space for M max = 1.8 Mq is marginal but increases rapidly as M max exceeds 2 Mq. 



iCook et al.lll994tlBaumgarte et al.ll2000f ) : the exact value 
of this limit depends on the radial profile of the differ- 
ential rotation. As a result, the delayed collapse channel 
requires that 

1.2 M max < Mi + M 2 < 1.4 Af max . (2) 

The second channel involves the prompt collapse of 
the remnant of the coalescence into a black hole and the 
formation of a > 0.1 Mq torus around it that is ac- 
creted onto the black hole in the required timescale. All 
numerical simulations of coalescing neutron stars show 
that equal mass mergers hamper the formation of these 
massive thick accretion disks. If, on the other hand, the 
ratio of the merging neutron star masses is larger than 
4:3, the lower mass neutron star fills its Roche lobe 
before contact and loses the required amount of mass 
into a torus before the remnant collapses into a black 
hole (jShibata fc Taniguchi|[2^MlDue9[2^Tol ). (The value 
of this ratio depends on the mass of the torus required 
to achieve the timescale of the GRBs by the subsequent 
accretion of this torus; see Rezzolla et al. 2010). This 
prompt collapse channel, therefore, requires that the to- 
tal mass of the coalescing neutron stars is > 1.4 M max 
and their mass ratio is > 4/3. 

In Figure 2, we calculate the outcome of mergers in 
the parameter space defined by the masses of the two 
coalescing neutron stars and delineate the areas in which 
short duration gamma-ray bursts can be generated via 
either channel. The two panels demonstrate the very 
strong sensitivity of these outcomes on the maximum 
mass that can be supported against collapse in a non- 
rotating neutron star. For a maximum neutron star mass 
of 1.8 Mq, the allowed area in this parameter space is 
marginal (left panel) but increases rapidly as M max grows 
to values larger than 2 Mq (right panel). Therefore, the 



unequivocal detection of a 1.97 Mq neutron star makes 
the idea of coalescing neutron stars as the central engines 
of short duration gamma-ray bursts viable. 

Figure 2 shows that the largest part of the param- 
eter space that leads to short duration gamma-ray 
bursts is populated by double neutron star systems in 
which one of the neutron stars is significantly more 
massive than the other. This is in contrast to the 
double neutron star systems observed in the Galactic 
disk, where both neutron stars have surp risingly simi- 
lar masses (jThorsett fc Chakrabartvl Il999l ). Indeed, a 
significant amount of mass transfer onto one of the two 
neutron stars is likely to take place only in a low-mass 
X-ray binary system, which does not lead to the forma- 
tion of a double neutron star. Double neutron stars with 
high mass ratios are perhaps most easily formed in glob- 
ular clusters via exchange interactions. Globular clusters 
have large populations of low-mass X-ray binaries and, 
therefore, may harbor the systems responsible for the 
short duration gamma-ray bursts (see also Grindlay et 
al. 2007). This would explain the recent observations of 
optical counterparts of short duration gamma-ray bursts, 
which appear to lie outside the host galaxies (e.g., Berger 
2010). 

4. MEASURING THE EQUATION OF STATE OF NEUTRON 
STARS WITH GRAVITATIONAL WAVES 

The coalescence of two neutron stars is also expected 
to be the primary source of gravitational waves de- 
tected with current and future ground- based detectors 
such as LIGO, VIRGO, and GEO 600 (|Anderson et al.l 
2009). Detailed modeling of the waveforms during the 
final stages of the coalescence can place strong con- 
stra ints on the equa tion of state of neutron-star mat- 
ter (|Read et al.ll2009T) . This will be possible with a gen- 
eration of detectors beyond Advanced LIGO, as it re- 
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quires following the gravitational wave signal to frequen- 
cies > 1 kHz. However, even the low-frequency wave- 
forms may encode the characteristics of the equation 
of state as they are affected by the tidal deformation 
experi enced by the two neutron stars at larger separa- 
tions ([Flanagan fc Hindered 120081 ) . 

The degree of tidal deformation, as measured by 
the Love number, is relatively insensitive to the com- 
pactness of the neutron star but depends strongly 
on how centr ally condensed t he density profile inside 
the stars is (IHinderer I 120081 : IDamour &; Nagarl 120091 : 
IBinnington fc Poissonll2009D . because stars that are cen- 
trally condensed are weakly deformed in the presence of 
an external tidal gravitational field. 

Equations of state that allow pulsar masses > 2Af Q 
typically yield neutron stars that are not centrally con- 
densed (as can be inferred from the polytropic mod- 
els of Hinderer 2008). Therefore, the discovery of the 



1.97 Mq pulsar indicates that coalescing neutron stars 
will undergo significant tidal deformation, and informa- 
tion about the equation of state of the stellar interior will 
be encoded in low- frequency (< 600 Hz) gravitational 
waves emitted during coalescence. Detection of these 
waves will allow accurate measurements of the equation 
of state of neutron-star matter in the near future. 
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